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The high spectral resolution and sensitivity of HerschelfHlFl allows for the detection of multiple rotational water lines and accurate determinations 
of water production rates in comets. In this letter we present HIFI observations of the fundamental lio-loi (557 GHz) ortho and In-Ooo (1113 
GHz) para rotational transitions of water in comet 8 IP/Wild 2 acquired in February 2010. We mapped the extent of the water line emission with 
five point scans. Line profiles are computed using excitation models which include excitation by collisions with electrons and neutrals and solar 
infrared radiation. We derive a mean water production rate of 1.0 x 10-** molecules s"' at a heliocentric distance of 1.61 AU about 20 days before 
perihelion, in agreement with production rates measured from the ground using observations of the 18-cm OH lines. Furthermore, we constrain 
the electron density profile and gas kinetic temperature, and estimate the coma expansion velocity by fitting the water line shapes. 
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1. Introduction 

Water is the main volatile constituent of cometary comae. Thus 
it is crucial to determine the water production rate to understand 
cometary activity and determine the relative abundances of other 
volatiles in the coma (see e.g. Bockelee-Morvan 1987). Water 
molecules in cometary atmospheres are excited by collisions 
with neutrals and electrons, and radiative pumping of the fun- 
damental vibrational levels by the solar infrared radiation. The 
fundamental lio-loi ortho-water transition at 556.936 GHz is 
one of the strongest lines in cometary spectra. This transition 
has been detected using the Submillimeter Wave Astronomical 
Satell ite (SWAS) (|Neufeld et al||2000|i, and subsequen tly with 



Odin ( ILecacheux et al.||2003j |Biver et al.|l2007[ [2009) 1 and the 
Herschel Space Observatory ( Hartogh et al.|2010| . Herschel pro- 
vided the first detection of the 2i2-loi (1669.905 GHz) ortho and 
lii-0()o (1113.343 GHz) para transitions of water in cometary 
atmospheres (Hart ogh et al.|2010 i. The continuum emission in 
comet C/2006 W3 (Christensen) has been measured by Herschel 
(Bockelee-Morvan et al. 201 ). 

Comet 8 IP/Wild 2 is a Jupiter family comet that was the tar- 



get of NASA's Stardust sample retuiTi mission (Brownlee et al. 
|2004| l. Thousands of particles in the submillimeter size range 
were collected for laboratory study during the Stardust flyby 
on January 2, 2004 ( ,Brownlee et al.n2006j . Optical images of 
the nucleus obtained by the navigation camera show very com- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



plex topographic features (Brownlee et al.'2004). Gas production 
rates for several molecules have been monitored at multiple ap- 
paritions (|Osip et al.|1992t|Makinen et al.|2001|[Fmk et al.|1999[ 



Farnham & Schleicher|2005) l. A maximum water production rate 
of 2h20 ~ 1-1 X 10^ s~' is estimated from narrowband pho- 
tometry of OH emission in the near-UV ( [Famham & Schleicher| 
[20051 ). 

The Heterodyne Instrument for the Far-Infrared (HIFI) on- 
board Herschel is well suited to observe water vapor in cometary 
coma with high spectral resolution (de Graauw et al. 2010} 



|Pilbrattetal.|2010| ). Comet 81P/Wild 2 was observed with HIFI 
on February 1^, 2010, in the framework of the Herschel guar- 
anteed ti me key project 'Water and related chemistry in the SoIm] 
Systemj ' ( [Hartogh et al.|2009| l. 



2. HIFI observations 

Comet 81P/Wild 2 passed its perihelion on February 22.7 UT at 
rji = 1 .598 1 AU from the Sun. The Herschel observations took 
place in early February about 20 days before perihelion, at he- 
liocentric distance rj, - 1.6 AU and a distance of A ^ 0.9 AU 
from Herschel. The comet could not be observed after February 
14 due to visibility constraints. HIFI can observe two polar- 
izations simultaneously. With a spectral resolution of 1.1 MHz 
(Wide Band Spectrometer; WBS) or 140 KHz (High Resolution 
Spectrometer; HRS) HIFI can resolve spectrally cometary line 
shapes and asymmetries ([ Hartogh et al.|2010l l. 

Table [T] shows a summary of the HIFI observations of the 
lio-loi (557 GHz, HIFI bands la and lb) and ln-Ooo (1113 
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GHz, HIFI band 4b) water lines in comet SlPAVild 2. The 
position of the comet and its relative motion with respect to 
Herschel were calculated using the JPL's HORIZONS system 
The Herschel telescope has a diameter of 3.5 m, with the cor- 
responding HIFI Half Power Beam Widths (HPBW) of 19.2" 
and 38.1" at 1113, and 557 GHz, respectively. The correspond- 
ing beam sizes projected on the comet are 1.5 and 2.2 x 10"* 
km, respectively. Observations were conducted in the standard 
dual beam switch (DBS) cross map observing mode with a chop- 
per speed of 4 Hz and a separation of the cross positions corre- 
sponding to a Nyquist sampling of the beam size (22", 19", and 
10" for observations in bands la, lb, and 4b, respectively). One 
axis of the cross was aligned along the Sun direction at PA - 
111°. The reference OFF position was 3' from the position of 
the comet. This observing mode was not released at the time of 
the observations. 



3. Data analysis 

3. 1 . Water line emission 

The data were reduced to the level 2 products using the stan- 
dard Herschel Interactive Processing Environment (HIPE) 3.0 
pipeline for HIFI ( |de Graauw et al.|2010[ ). Integrated line inten- 
sities and velocity offsets in the comet rest frame are given in 
Tablefllfor the central point in the maps. Line intensities are cal- 
culated from the weighted averages of spectra in horizontal (H) 
and vertical (V) polarizations. We scaled the main beam bright- 
ness temperature using the beam efficiency of the Herschel tele- 
scope of 0.7 and 0.75 for the 1113 and 557 GHz lines respec- 
tively (see e.g. Hartogh et al.|2010 l. The uncertainties in line in- 
tensity, velocity shift and water production rates are 1-cr statisti- 
cal uncertainties. Uncertainties in the sideband gain ratio, which 
reach 5% ( jRoelfsema et al.||2010| l, are considered in the data 
analysis. 

Figurefllshows the HRS spectrum of the lio-loi ortho-water 
line at 557 GHz observed on February 4.24 UT at the center of 
the map. The central spectrum of the map of the Ih-Oqo para- 
water line at 1113 GHz, averaging the two observations to in- 
crease the signal-to-noise ratio, is shown in Fig.|2] These spec- 
tra are uniformly weighted averages of the H and V polariza- 
tions. They correspond also to the average of the central points 
of the horizontal and vertical drifts of the cross map. The fre- 
quency scale of the observed spectra was corrected for the geo- 
centric velocity of the comet and the spacecraft orbital velocity. 
The lines are optically thick and slightly asymmetric due to self- 
absorption effects in the foreground coma. We estimate a mean 
expansion velocity Vexp = 0. 6 kms"' from the width of the 557 
GHz line ( jBiver et al.p007| . 



Expansion velocities in the range 
- 0.5 - 0.8 kms~^ are typical for weak comets. 
Mapping observations at 557 and 1113 GHz are shown in 
Figs. [3] and [4] and online Fig. [5] 



3.2. Radiative transfer modelling 

Molecular excitation in the outer coma is dominated by colli- 
sions with electrons and infrared fluorescence by solar radiation. 
We use a spherically symmetric Monte Carlo radiative transfer 
numerical code to compute the populations of the water rota- 
tional levels as function of the distance from the nucleus (see 
[Hartogh et al.|2010 and references therein). A constant gas tem- 
perature of r = 40 K is assumed. The ortho-to-para water ratio 
in comets varies from 2.5 to 3 (e.g. Bonev et al.||2007 and ref- 
erences therein). We assume a value of 3. The scaling factor of 
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Fig. 1. Central point in the DBS map of the ortho- water line 
at 556.936 GHz in comet 81P/Wild 2 obtained by the HRS on 
February 4.24 UT. The velocity scale is given with respect to the 
comet rest frame with a resolution of ~ 32 m s"' . A synthetic line 
profile for isotropic outgassing with Vexp = 0.6 km s"', T = 40 K, 
and x„^ - 0.2 is shown by the dashed line. 




Fig. 2. Average of the HRS observations of the para-water line 
at 1113.343 GHz towards the nucleus of comet 8 IP/Wild 2 ob- 
tained on February 2.13 UT and 2.15 UT. The velocity resolution 
is ~ 32 m s"' after smoothing. The dashed line shows the best fit 
profile for the same parameters as in Fig. [T] 



the electron density profile with respect to a nominal profile de- 
duced from IP/Halley in situ data is chosen to be x„^ = 0.2 and 
1 . The former value was found to best explain the brightness dis- 
tribution of the 1 lo-loi line obtained from mapping observations 
( |Biver et al.|2"007| [Hartogh et al.|2010| l. 

The radiative transfer equation is solved along different lines 
of sight through the coma covering 2.5 x HPBW, and the beam 
averaged emission spectra are computed. We obtain the water 
production rate by comparing the calculated and observed line 
intensities. The synthetic line profiles agree approximately with 
the observed line shapes (Figs. [l]-^. The differences (overall 
shift and excess emission in the blue wing of the profiles) suggest 
preferential outgassing towards the Sun and a day/night asym- 
metry in the gas velocity field. 
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Table 1. HIFI observations of water vapour in comet 8 1 P/Wild 2 and derived water production rates. The observations were acquired 
in February 2010 using the standard fast DBS cross map observing mode. The intensity, velocity shift and QhiO correspond to the 
average of the central points in the cross maps. 



Obs. ID 


UT start date 


A 


rh 


Band 


Frequency 


Integration time 


Intensity 


Velocity shift 


eH,o" 




[mm/dd.ddd] 


[AU] 


[AU] 




[GHz] 


[s] 


[Kkms-'] 


[ms-'] 


[10^8 s-'] 


1342190098 


02/1.443 


0.949 


1.612 


la 


556.9360 


684 


1.807 ±0.012 


4±5 


1.00 ±0.05 


1342190185 


02/4.235 


0.927 


1.609 


lb 


556.9360 


748 


2.077 ± 0.015 


-10±5 


1.13 ±0.06 


1342190231 


02/2.134 


0.943 


1.611 


4b 


1113.3430 


390 


1.997 ±0.053 


-58±21 


1.07 ±0.06 


1342190232 


02/2.151 


0.943 


1.611 


4b 


1113.3430 


366 


1.594 ±0.052 


-41 ±26 


0.86 ± 0.05 



(n) Water production rates computed with T = 40 K and x„^ = 0.2. The error bar includes a 5% calibration uncertainty in the sideband gain ratio. 
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Fig. 3. DBS cross-map of the ortho-water lio-loi transition ob- 
tained with the HRS on February 1 .44 UT. Offset positions are 
~ 22" apart. The velocity resolution is ~ 32 ms"'. The Sun 
direction is indicated. 
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Fig. 4. Average of the cross-maps of the 1 1 i-Ooo para-water line 
obtained with the HRS on February 2.13 and 2.15 UT. The off- 
set spacing is ~ 10" and the effective velocity resolution after 
smoothing is ~ 32 ms"'. The Sun direction is indicated. 



3.3. Water production rates 

Mapping observations can be used to constrain excitation pa- 
rameters such as the neutral gas kinetic temperature and electron 
density scaling factor The spacing of 10" and 22" between the 
observed positions corresponds to distances of ~ 1.2 x 10"* and 
2.4 X lO'* km at the distance of the comet, respectively, where 
the excitation of the water rotational levels is controlled both by 
colUsions with electrons and infrared fluorescence (see Fig. 5 in 
[Hartogh et al.|20T0) l. Figure [6] shows the evolution of the line 
area of the 557 and 1113 GHz lines as a function of the posi- 
tion offset. As shown in Fig. |6j model calculations performed 
with r = 40 K and x„^ = 0.2 provide a satisfactory fit to the 
observed radial brightness profiles. For a larger electron content 
(x„^ -I), the line intensities at the offset positions are lower 
than observed. Production rates for x„^ = 0.2 are consistent with 
values obtained from OH observations (Crovisier et al. in prepa- 
ration), while those for x„ . = 1 are about a factor of two lower 

The water production rates derived from the 557 and 1113 
GHz lines are consistent with each other and in the range 
0.9-1.1 X 10^^ s"'. IVIodel calculations show that the 7(557 
GHz)//(1113 GHz) intensity ratio does not depend strongly on 



the x„^ parameter. Gas kinetic temperatures T in the range 30- 
50 K provide a satisfactory fit of the measured intensity ratio, 
averaging the lines observed at the different dates and consider- 
ing uncertainties in flux calibration. Lower temperatures (by ~ 
15 K) are retrieved when the ortho-to-para ratio is taken equal 
to 2.5. Such low T values are typical for comets of relatively 
low activity. However, these temperature retrievals are possibly 
affected by the non-simultaneity of the observations. 

3.4. Comparison witli observations of water 
pliotodissociation products 

In support of the Herschel observations, the OH lines at 18-cm 
were observed with the Nangay radio telescope from November 
16, 2009 to February 14, 2010 (Crovisier et al. in prepara- 
tion). The methods of observation and analysis are described in 
ICrovisier et al. ( 2002| l. The comet could be detected only after 
integrating over many days in December and January. It could 
not be detected close to perihelion because the OH maser inver- 
sion was small at that time. Water is the only significant par- 
ent molecule of the OH radical in cometary atmospheres and 
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Fig. 6. Intensities of the lio-loi (Feb. 1.44 UT, left panel) and 
lii-Ooo (Feb. 2.13-2.15 UT, right panel) lines as a function of 
position offset from the estimated position of the peak bright- 
ness. Model calculations for x„. = 0.2 and 1 are shown with a 
black solid line and red dashed line, respectively. We assumed T 
= 40 K and Vexp = 0.6 km s"' . 



the lifetimes of the two species are well known, allowing for a 
determination of the water production rate (see e.g. Cochran & 
|Schleicher| 1 993| l. The retrieved water production rates, includ- 
ing the HIFI retrievals, are shown in Fig. [7] 

Figure |7] includes water production rates measured for the 
1997 apparitions from OH narrowband photometry ( Farnham & 
Schleicher]|2005 i, spectroscopy of the 'D line of atomic oxy- 



gen ( Fink et al.||1999| l and Lyman-a measurements made with 
the SWAN instrument onboard SOHO (Miikinen et al. 2001). 
The latter water production rates are about a factor of two higher 
than those obtained from narrowband photometry of OH emis- 
sion ( Farnham & Schleicher|2003 1. The higher production rates 
measured 80-20 days before perihelion may indicate a seasonal 
effect, with a water production rate peak at ~ 40-60 days pre- 
perihelion (see e.g. Szutowicz et al.|2()08 1. 



4. Conclusions 

Although comet 8 IP/Wild 2 has been extensively observed from 
the ground, it is still poorly observed at sub-millimeter wave- 
lengths. 81P/Wild 2 was observed with HIFI in the period 1-4 
February 2010 at 1.6 AU heliocentric distance and 0.95 AU from 
Herschel. The fundamental ortho- and para-water rotational 
transitions at 556.936 GHz and 1113.343 GHz were detected 
at high spectral resolution. DBS cross-maps were obtained to 
study the excitation conditions throughout the coma. We calcu- 
late a water production rate in the rangeO. 9-1. 1x10^*^ s"' using a 
radiative transfer code which includes collisional effects and in- 
frared fluorescence by solar radiation. Water production rates are 
in good agreement with those derived from ground-based obser- 
vations of the OH 18-cm emission at the Nanfay radio telescope 
about 30 days before perihelion. 
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Fig. 7. Water production rates with 1-cr uncertainties in comet 
8 IP/Wild 2 as a function of time from perihelion. HIFI measure- 
ments are shown by squares. Filled circles represent OH 18-cm 
line observations with the Nan^ay radio telescope. Water pro- 
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Fig. 5. DBS cross-map of the lio-loi ortho-water Hne obtained 
with the HRS on February 4.24 UT. The offset spacing is ~ 19" 
and the velocity resolution is ~ 32 ms"'. The Sun direction is 
indicated. 



